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ABSTRACT 

We present detailed imaging of Faraday rotation and depolarization for the radio galaxies 
0206+35, 3C 270, 3C 353 and M 84, based on Very Large Array observations at multiple fre- 
quencies in the range 1365 to 8440 MHz. All of the sources show highly anisotropic banded 
rotation measure (RM) structures with contours of constant RM perpendicular to the major 
axes of their radio lobes. All except M 84 also have regions in which the RM fluctuations have 
lower amplitude and appear isotropic. We give a comprehensive description of the banded 
RM phenomenon and present an initial attempt to interpret it as a consequence of interactions 
between the sources and their surroundings. We show that the material responsible for the 
Faraday rotation is in front of the radio emission and that the bands are likely to be caused 
by magnetized plasma which has been compressed by the expanding radio lobes. We present 
a simple model for the compression of a uniformly-magnetized external medium and show 
that RM bands of approximately the right amplitude can be produced, but for only for special 
initial conditions. A two-dimensional magnetic structure in which the field lines are a family 
of ellipses draped around the leading edge of the lobe can produce RM bands in the correct 
orientation for any source orientation. We also report the first detections of rims of high de- 
polarization at the edges of the inner radio lobes of M 84 and 3C 270. These are spatially 
coincident with shells of enhanced X-ray surface brightness, in which both the field strength 
and the thermal gas density are likely to be increased by compression. The fields must be 
tangled on small scales. 

Key words: - galaxies: magnetic fields - radio continuum: galaxies - (galaxies:) intergalactic 
medium - X-rays: galaxies: clusters 



1 INTRODUCTION 

The detection of diffuse synchrotron emission (radio halos) on Mpc 
scales in an increasing number of galaxy clusters provides good 
evidence for a distributed magnetic field of /xGauss st rength in 
the hot intracluster medium (ICM; see e.g. Fer rari et al. 2008 for 
a review). Imaging of Faraday rotation of linearly-polarized radio 
emission from embedded and background sources confirms that 
there are fields associated with thermal plasma a long lines of sight 
through the clusters (e.g. ICarilli & Tavloil 12002). Observations of 
Faraday rotation can also be made for radio galaxies in sparser en- 
vironments, allowing the study of magnetic field s in environments 
too sparse for radio halos to be detected (e.g. lLaing et al.l [20081 : 

iGuidetti et ai1l2010h . 

The Faraday effect (Faraday] 1 18461) is the rotation suffered 



* E-mail: dguidett@eso.org 



by linearly polarized radiation travelling through a magnetized 
medium, and can be described by the two following relations: 

A*[ rad ] = *(A)[ rad ] - * [ rad j = A 2 m 2] RM [radm -2], (1) 

with 

RM[r adm -2] = 812 / n c [cm -3jB 2 [ MG ]dz[ kpc ] , (2) 
Jo 

where ^(A) and are the E-vector position angle of linearly po- 
larized radiation observed at wavelength A and the intrinsic angle, 
respectively, n e is the electron gas density, B z is the magnetic field 
along the line-of-sight (-By), and L is the integration path. RM is 
the rotation measure. 

Observations of Faraday rotation variations across extended 
radio galaxies allow us to derive information about the integral of 
the density-weighted line-of-sight field component. The hot (T ~ 
10 7 — 10 8 K) plasma emits in the X-ray energy band via thermal 
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bremsstrahlung. When high quality X-ray data for a radio-source 
environment is available, it is possible to infer the gas density dis- 
tribution and therefore to separate it from that of the magnetic field, 
subject to some assumptions about the relation of field strength and 
density. 

Most of the RM images of radio galaxies published so far 
show patchy structures with no clear preferred direction, consis- 
tent with isotropic foreground fluctuations over a r ange of linear 
scales ranging from tens of kpc to < 100 pc (e.g. Govoni et al. 
l2006l: IGuidetti et al1l2008l : lLaing et aljl2008l ; IGuidetti et alj|2010h . 
Numerical modelling has demonstrated that this type of complex 
RM structure can be accurately reproduced if the magnetic field 
is randomly variable with fluctuations on a wide range of spatial 
scales, and i s spread throughout the whole group or cluster envi- 



ronm e nt (e.g. Murgia et alj2004 ; Govoni et al.l2006l:lGuidetti et al.l 
120081 : lLaing et al] 120081 ; IGuidetti et alj l20ld : IVacca et al.l l2010h . 
These authors used forward modelling, together with estimators of 
the spatial statistics of the RM distributions (structure and auto- 
correlation functions or a multi-scale statistic) to estimate the field 
strength, its relation to the gas density and its power spectrum. The 
technique of Bayesian maximum likelihood has also been u sed for 
this purpose fenfilin & Vogtll2005l;lKu~char & En61inll2009l) . 

In order to derive the three-dimensional magnetic field power 
spectrum, all of these authors had to assume statistical isotropy for 
the field, since only the component of the magnetic field along the 
line-of-sight contributes to the observed RM. This assumption is 
consistent with the absence of a preferred direction in most of the 
RM images. 

In contrast, the present paper reports on anisotropic RM struc- 
tures, observed in lobed radio galaxies located in different environ- 
ments, ranging from a small group to one of the richest clusters of 
galaxies. The RM images of radio galaxies presented in this paper 
show clearly anisotropic "banded" patterns over part or all of their 
areas. In some sources, these banded patterns coexist with regions 
of isotropic random variations. The magnetic field responsible for 
these RM patterns must, therefore, have a preferred direction. 

One source whos e RM structure is dom inated by bands is al- 
ready known: M 84 jLaing&Bridie1ll987h . In addition, there is 
some evidence for RM bands in source s which also show stro ng 
irregular fluctuations, such as Cygnus A dCarilli & Tavlonl2002l) . It 
is possible, however, that some of the claimed bands could be due 
to imperfect sampling of an isotropic RM distribution with large- 
scale power, and we return to this question in Section [9~Tl 

The present paper presents new RM images of three sources 
which show spectacular banded structures, together with improved 
data for M 84. The environments of all four sources are well char- 
acterized by modem X-ray observations, and we give the first com- 
prehensive description of the banded RM phenomenon. We present 
an initial attempt to interpret the phenomenon as a consequence of 
source-environment interactions and to understand the difference 
between it and the more usual irregular RM structure. 

The RM images reported in this paper are derived from new 
or previously unpublished archive Very Large Array (VLA fl data 
for the nearby radio galaxies 0206+35 , M84 fcaing et all 1201 ll) . 
3C 270 jLaing, Guidetti & Bridld20 111) and 3C 353 (Swain, private 
communication: see lSwainll 19961) . 

The paper is organized as follows. In Section© the radio and 
X-ray properties of the sources under investigation are presented 
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and in Section[3]we briefly summarize the techniques used to anal- 
yse depolarization and two-dimensional variations of RM. Sec- 
tions!?] and [5] present the RM and depolarization images on which 
our analysis is based and correlations between the two quantities. In 
Section[6] we evaluate the RM structure functions in regions where 
the fluctuations appear to be isotropic and derive the power spec- 
tra. A simple model of the source-environment interaction which 
characterises the effects of compression of a magnetised IGM is 
described in Section [7] This can produce RM bands, but only un- 
der implausible special initial conditions. Empirical "draped" field 
configurations which are able to reproduce the banded RM distri- 
butions are investigated in Section[8] In Section[9] we speculate on 
correlations between radio source morphology, and RM anisotropy, 
discuss other examples from the literature, consider the effects of 
an isotropic foreground Faraday screen on the detectability of RM 
bands and briefly discuss possible asymmetries in the amplitude of 
the RM bands between the approaching and receding lobes. Finally, 
SectionfTOl summarizes our conclusions. 

Throughout this paper we assume a ACDM cosmology with 
Ho = 71 km s^Mpc" 1 , « ro = 0.3, and fi A = 0.7. 



2 THE SAMPLE 

High quality radio and X-ray data are available for all of the 
sources. In this Section we summarize those of their observational 
properties which are relevant to our RM study. A list of the sources 
and their general parameters is given in Table[T] while Table© 
shows the X-ray parameters taken from the literature and equipar- 
tition parameters derived from our radio observations. 

The sources were observed with the VLA at several frequen- 
cies, in full polarization mode and with multiple configurations so 
that the radio structure is well sampled. The VLA observations, 
data reduction and detailed descr iptions of the rad io structures are 
given for 0206+35 and M 84 by lLaing et all d201ll) for 3C 270 b y 
lLaing, Guidetti & Bridle! ( 1201 if) , and for 3C 353 bv lSwairj dl996h . 
All of the radio maps show a core, two sided jets and a double- 
lobed structure with sharp brightness gradients at the leading edges 
of both lobes. The synchrotron minimum pressures are all signif- 
icantly lower than the thermal pressures of the external medium 
(Table©. 

All of the sources have been observed in the soft X-ray band 
by more than one satellite, allowing the detection of multiple com- 
ponents on cluster/group and sub-galactic scales. The X-ray mor- 
phologies are characterized by a compact source surrounded by ex- 
tended emission with low surface brightness. The former includes a 
non-thermal contribution, from the core and the inner regions of the 
radio jets and, in the case of 0206+35 and 3C 270, a thermal com- 
ponent which is well fitted by a small core radius (3 model. The 
latter component is associated with the diffuse intra-group or intra- 
cluster medium. Parameters for all of the thermal components, de- 
rived from X-ray observations, are listed in Table© Because of the 
irregular morphology of the hot gas surrounding 3C 353 and M 84, 
it has not been possible to fit /3 models to their X-ray radial surface 
brightness profiles. 



2.1 0206+35 (4C 35.03) 

0206+35 is an exte nded Fanaroff-Riley Class I (FRI; 
iFanaroff & Rileyl 1 19741) radio source whose optical counter- 
part, UGC1 1651, is a D-galaxy, a member of a dumb-bell system 
at the centre of a group of galaxies. At a resolution of 1.2 arcsec the 
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Figure 1. X-ray images overlaid with radio contours for all sources: (a) 0206+ 35: 1385.1 MHz VLA radio m ap with 1.2 arcsec FWHM; the contours are spaced 
by factor of 2 between 0.06 and 15 mjy beam" 1 . The ROSAT PSPC image dWorrall & Birkins haw 2000), is smoothed with a Gaussian of a = 30 arcsec. 
(b) M 84: 1413.0 MHz VLA radio map with 4.5 arcsec FWHM; the contours are spaced by factor of 2 between 1 and 128 mjy beam -1 . The Chandra 
iFinoguenov & Jones 11200 lb image is a wavelet reconstruction on angular scales from 4 up to 32 arcsec. (c) 3C270: 1365.0 MHz VLA radio map with 
5.0 ar csec FWHM; the contours are spaced by factor of 2 between 0.45 and 58 mjy beam -1 . The XMM-Newton image )Croston."H ardcastle & Birki nshaw I 
|2005|) is smoothed with a Gaussian of a = 26 arcsec. (d) 3C 353: 1385.0MHz VLA radio map with 1.3 arcsec FWHM; the contours are spaced by factor of 
2 between 0.35 and 22 mjy beam -1 . The XMM-Newton image t Goodger et al. 2008?) is smoothed with a Gaussian of a = 30 arcsec. The Chandra image of 
3C 270 is displayed in logarithmic scale. 



Table 1. General optical and radio properties: Col. 1: source name; Cols. 2&3: position; Col. 4: redshift; Col. 5: conversion from angular to spatial scale with 
the adopted cosmology; Col. 6: Fanaroff-Riley class; Col. 7: the largest angular size of the radio source; Col. 8: radio power at 1.4 GHz; Col. 9: angle to the 
line of sight of the jet axis; Col. 10: radio spectral index; Col. 11: environment of the galaxy; Col. 12: reference. 



source 


RA 


DEC 


z 


kpc/arcsec 


FR class 


LAS 


log Pi A 


e 


env. 


ref. 




[J2000] 


[J2000] 








[arcsec] 


[W Hz" 1 ] 


[degree] 






0206+35 (4C 35.03) 


02 09 38.6 


+35 47 50 


0.0377 


0.739 


I 


90 


24.8 


40 


group 


1 


3C353 


17 20 29.1 


-00 58 47 


0.0304 


0.601 


II 


186 


26.3 


90 


poor cluster 


2 


3C270 


12 19 23.2 


+05 49 31 


0.0075 


0.151 


I 


580 


24.4 


90 


group 


3 


M84 


12 25 03.7 


+12 53 13 


0.0036 


0.072 


I 


150 


23.2 


60 


rich cluster 


3 



References for the environmental classification: (1) Miller et al. '2002); (2) de Vaucoleurs i 1991): (3) Traser et al '2000). 
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Table 2. X-ray and radio equipartition parameters for all the sources. Col. 1: source name; Col. 2: X-ray energy band; Col. 3: average thermal temperature; 
Cols. 4,5,6 and 7,8,9 best-fitting core radii, central densities and j3 parameters for the outer and inner /3 models, respectively; Col. 10: average thermal pressure 
at the midpoint of the radio lobes; Cols. 1 1&12: minimum synchrotron pressure and corresponding magnetic field; Col. 13: references for the X-ray models. 



source band kT r CXout no out Pout r cx in n 0in /3 in P P m in -Bpmin ref. 

[keV] [keV] [kpc] [cm -3 ] [kpc] [cm" 3 ] [dyne cm -2 ] [dyne cm -2 ] /iG 



0206+35 


0.2-2.5 


1 3+ 1 ' 3 
1 "'-0.3 


22.2 


2.4 x 


10- 3 


0.35 


0.85 


0.42 


0.70 


9.6 x 10" 


12 


4.31 x itr 


13 


5.70 


1, 2 


3C353 




4 33+0-25 
*- 33 -0 .24 




















1.66 x 10" 


-12 


11.2 


3 


3C270 


0.3-7.0 


1 45+O.23 


36.8 


7.7 X 


10-3 


0.30 


1.1 


0.34 


0.64 


5.75 x 10- 


-12 


1.64 x 10- 


-13 


3.71 


4 


M84 


0.6-7.0 


Q6 +0.05 
u -°-0.05 










5.28±0.08 


0.42 


1.40±0.03 


1.70x 10" 


-11 


1.07x 10" 


-12 


9.00 


5 



References: (I) Wonall & Birkinshaw i 2000); (2 iWonall, Birkinshaw & Hardcastle i 2001); (3) Iwasawa et al. 2000); (4) Croston et al. i 2008); (5) i Finoeuenov & Jones 2001). 



radio emission shows a core, with smooth two-sided jets aligned in 
the N W-SE direction a nd surr ounded by a diffuse and symmetric 
halo. lLaing & Bridle] fcOllbl) have estimated that the jets are 
inclined by ~40° with respect to the line of sight, with the main 
(approaching) jet in the NW direction. 

0206+35 has been observed wit h both the ROSAT 
PSPC and HRI instruments jWorrall & Birkinshaw! 
1 1994 l200(jl: iTrussoni et al.1 Il997h and with Chandra 
dWorrall. Birkinshaw & Hardcastld 1200 lh . The X-ray emission 
consists of a compact source surrounded by a galactic atmosphere 
which merges into the much more extended intra-group gas. The 
radius of the extended halo observed by the ROSAT PSPC is 
?s2.5 arcmin (Fig. [TJi). The ROSAT and Chandra X-ray surface 
brightness profiles are well fit by the combination of f} models with 
two different core radii and a power-law component (Hardcastle, 
private communication; Table©. 



mate for the inclination of the jets is ~90° ( Swain, Bridle & Bauml 
1998). The eastern jet is slightly brighter and ends in a well-defined 
hot spot. The radio lobes have markedly different morphologies: 
the eastern lobe is round with sharp edges, while the western lobe is 
elongated with an irregular shape. The location of the source within 
the cluster is of particular interest for this work and might account 
for the different shapes of the lobes. FigJjJd) shows the XMM- 
Newton image overlaid on the radio contours. The image shows 
only the NW part of the cluster, but it is clear that the radio source 
lies on the edge of the X-ray emitting gas distribution, so that the 
round eastern lobe is encountering a higher external density and is 
probably also behind a larger column of F araday-rotating material 
dlwasawa~et al. 2000, Goodger et alj2008h . In particular, the image 
published bv lGoodger et al. (2008) shows that the gas density gra- 
dient persists on larger scales. 



2.2 3C270 

3C 270 is a radio source classified as FRI in most of the literature, 
although in fac t, the two lobes have differen t FR classifications at 
low resolution ( lLaing. Guidetti & Bridlell201 J). The optical coun- 
terpart is the giant elliptical galaxy NGC 4261, located at the centre 
of a nearby group. The radio source has a symmetrical structure 
with a bright core and twin jets, extending E-W and completely 
surrounded by lobes. The low jet/counter-jet ratio indicates that the 
jets are cl ose to the plane of the sky, wi th the Western side ap- 
proaching ( lLaing. Guidetti & BridleEoilh . 

The XMM-Newton image (Fig.QJ) shows a disturbed dis- 
tribution with regions of low surface-brightness (cavities) at 
the positi ons of both radio lobes. A recent Chandra ob- 
servation dWorrall et al] l2010h shows "wedges" of low X- 
ray surface brightness surrounding the inner jets (see also 
Crosto n. Hardcastle & Birkinshaw 2005, Finog uenov et alj|2006l , 
Jetha et alj|2007l ICroston et alj|2008h . The overall surface bright- 
ness profile is accurately reproduced by a point source con- 
volved wi th the Chandra p oint spread func tion plus a double 
/3 model dCroston et al] 120081 projb model). ICroston et aDl2008l 
found no evidence for a temperature gradient in the hot gas. The 
group is characterized by high temperature and low luminosity 
dFinoguenov et al.l2006T) . which taken together provide a very high 
level of entropy. This might be a further sign of a large degree of 
impact of the AGN on the environment. 

2.3 3C353 

3C 353 is an extended FRII radio source identified with a D-galaxy 
embedded at the periphery of a cluster of galaxies. The best esti- 



2.4 M84 

M84 is a giant elliptical galaxy located in the Virgo Cluster at 
about 400 kpc from the core. Optical emission-line imaging shows 
a disk of ionized gas aroun d the nucleus, with a maximum detected 
extent of 20 x 7arcsec 2 {Hansen, Norgaard-Nielsen & Jorgensenl 
ll985l : lBaum et alJl988l ; lBower et al.ll 1997ll2000h . The radio emis- 
sion of M 84 (3C 272.1) has an angular extension of about 3 arcmin 
(~ 1 1 kpc) and shows an unresolved core in the nucleus of the 
galaxy, two resolved jet s and a pair of wide lobes (Laing & Bridle 
1 19871 : lLaing et al-llzOlTh . The inclination to the line-of-sight of the 
inner jet axis is ~60° , with the northern jet approaching, but there 
is a noticeable bend in the coun ter-jet very close to t he nucleus, 
which complicates modelling (Laing & Bridle 2 01 lbl) . After this 
bend, the jets remain straight for ~40 arcsec, then both of them 
bend eastwards by ~90° and fade into the radio emission of the 
lobes. 

The morphology of the X-ray emission has a H-shape made up 
of shells of compress ed gas surrounding cavities coincident with 
both t he radio lobes (Finoguen ov et al. I [2008; Finoguen ov et al] 
2006; iFinoguenov & Jones I 1200 lh . This shape, together with the 
fact that the initial bending of the radio jets has the same direc- 
tion and is quite symmetrical, suggests a combination of interaction 
with the radio plasma and m otion of the galaxy within the cluster 
(Finogueno v & Jones 1 12001). The ratio between the X-ray surface 
brightness of the shells of the compressed gas and their surround- 
ings is ~3 and is almost constant around the source. The shells are 
regions of enhanced pressure and density and low entropy: the am- 
plitude of the density enhancements (a factor of ~3) suggests that 
they are produced by weak sho ck waves (Mach numbe r M ~ 1.3) 
driven by the expanding lobes ^Finoguenov et all2 006). 
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Figure 2. RM images for all sources: (a) 0206+35; (b) M 84; (c) 3C 270; (d) 3C 353. The angular resolution and the linear scale of each map are shown in the 
individual panels. 



3 ANALYSIS OF RM AND DEPOLARIZATION IMAGES 

For a fully-resolved foreground Faraday screen, the A 2 relation of 
Eq.[2]holds exactly and there is no change of degree of polarization, 
p, with wavelength. Even in the presence of a small gradient of 
RM across the beam, A 2 rotation is observed over a wide range of 
polarization angle. In this case, the emission tends to depolari ze 
with increasing wavelength, following the Burn law teurrJll966h : 

p(A) =p(0)exp(-fcA 4 ), (3) 

where p(0) is the intrinsic value of the degree of polarization 
and k=2\VRM\ 2 a 2 , with FW HM = 2<j(21n2) 1/2 iBurnll 19661 : 
lTribbldll99ll : lLaing et al.ll2008h - 



Since k oc \VRM\ 2 , Eq.[3] clearly illustrates that higher RM gra- 
dients across the beam generate higher k values and in turn higher 
depolarization. 

Our observational analysis is based on the following proce- 
dure. We first produced RM and Burn law k images at two differ- 
ent angular resolutions for each source and searched for regions 
with high k or correlated RM and k values, which could indicate 
the presence of internal Faraday rotation and/or strong RM gradi- 
ents across the beam. In regions with low k where the variations of 
RM are plausibly isotropic and random, we then used the structure 
function (defined in Eq.[4j to derive the power spectrum of the RM 
fluctuations. Finally, to investigate the depolarization in the areas of 
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Figure 3. Burn law k images for all sources at the same angular resolutions as for the RM images: (a) 0206+35 at 1 .2 arcsec FWHM; (b) M 84 at 4.5 arcsec 
FWHM; (c) 3C 270 at 5 arcsec FWHM; (d) 3C 353 at 1 .3 arcsec FWHM. The corresponding integrated depolarization (DP; Section[5} is indicated on the top 
right angle of each panel. The colour scale is the same for all displays. 



isotropic RM, and hence the magnetic field power on small scales, 
we made numerical simulations of the Burn law k using the model 
power spectrum with different minimum scales and compared the 
results with the data. 

The structure function is defined by 



S(r ± ) = < [RM(rj_ + r' x ) - RM(r x )] 2 > 



(4) 



dSimonetti. Cordes & Spangierl 1 19841 : iMinter & Spangled 1 19961) 
where r± and r' ± are vectors in the plane of the sky and () is an 
average over r j_ . 



We assume RM power spectra of the form: 

C(f±) = CofZ q /x</max 

= /_!_ > /max 

where f± is a scalar spatial frequency and fit the observed struc- 
ture function (including the effect of the observing beam) u sing the 
Hankel-transform method described by lLaing et alj (]2008) to de- 
rive the amplitude, Co and the slope, q. To constrain the RM struc- 
ture on scales smaller than the beamwidth, we estimated the mini- 
mum scale of the best fitted field power spectrum, A m in = 1/ /max, 
which predicts a mean value of k consistent with the observed one. 
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Table 3. Frequencies, bandwidths and angular resolutions used in the RM 
and Bum law k images discussed in Sects. |4]and|5] respectively. 



source 


V 


Au 


beam 




[MHz] 


[MHz] 


[arcsec] 


Ui-UUTJJ 


1385 1 


25 


1 2 




1464 9 


— j 






4885 1 


50 






1385 


12 5 


1,3 




1665 


12 5 






4866 3 


1 9 5 






8439.9 


12.5 




3C270 


1365.0 


25 


1.65 




1412.0 


12.5 






4860.1 


100 






1365.0 


25 


5.0 




1412.0 


12.5 






1646.0 


25 






4860.1 


100 




M84 


1413.0 


25 


1.65 




4885.1 


50 






1385.1 


50 


4.5 




1413.0 


25 






1464.9 


50 






4885.1 


50 





In this paper, we are primarily interested in estimating the RM 
power spectrum over limited areas, and we made no attempt to de- 
termine the outer scale of fluctuations. 

The use of the structure function together with the Burn law 
k represents a powerful technique to investig ate the RM power 
spectrum over a w ide range of spatial scales iLaing etalj 2008; 
lGuidettietal.l2010h . The two quantities are complementary, in that 
the structure function allows us to determine the power spectrum of 
the fluctuations on scales larger than the beamwidth, while the Burn 
law k constrains fluctuations of RM below the resolution limit. 



4 ROTATION MEASURE IMAGES 

The RM images and associated rms errors were produced by 
weighted least-squares fitting to the observed polarization angles 
^(\) as a fu nction of A 2 (Eq. Tj a t three or four frequencies (Ta - 
ble[3] see also lLaing et alj201l[ and lLaing, Guidetti & Bridld201lh 
using the RM task in the AIPS package. 

Each RM map was calculated only at pixels with rms 
polarization-angle uncertainties <10° at all frequencies. We re- 
fer only to the lower-resolution RM and k images for 3C 270 and 
M84 (Table [3}, as they show more of the faint, extended regions 
of these sources and are fully consistent with the higher-resolution 
versions. The RM ima ge of M 84 is consistent with that shown by 
lLaing&BridielJl987h , but is derived from four-frequency data and 
has a higher signal-to-noise ratio. 

In the fainter regions of 0206+35 (for which only three fre- 
quencies are available and the signal-to-noise ratio is relatively 
low), the RM task occasionally failed to determine the n7r ambi- 
guities in position angle correctly. In order to remove these anoma- 
lies, we first produced a lower-resolution, but high signal-to-noise 
RM image by convolving the 1 .2 arcsec RM map to a beamwidth of 

5 arcsec FWHM. From this map we derived the polarization-angle 
rotations at each of the three frequencies and subtracted them from 
the observed 1.2 arcsec polarization angle maps at the same fre- 



quency to derive the residuals at high resolution. Then, we fit the 
residuals without allowing any n7r ambiguities and added the result- 
ing RM's to the values determined at low resolution. This procedure 
allowed us to obtain an RM map of 0206+35 free of significant de- 
viations from A 2 rotation and fully consistent with the 1.2-arcsec 
measurements. 

We have verified that the polarization angles accurately follow 
the relation A^f oc A 2 over the full range of position angle essen- 
tially everywhere except for small areas around the optically-thick 
cores: representative plots of $ against A 2 for 0206+35 are shown 
in Fig.|4] The lack of deviations from A 2 rotation in all of the radio 
galaxies is fully consistent with our assumption that the Faraday 
rotating medium is mostly external to the sources. 

The RM maps are shown in Fig.|2] The typical rms error on 
the fit is ~2radm~ 2 . No correction for the Galactic contribution 
has been applied. 

All of the RM maps show two-dimensional patterns, RM 
bands, across the lobes with characteristic widths ranging from 3 
to 12 kpc. Multiple bands parallel to each other are observed in the 
western lobe of 0206+35, the eastern lobe of 3C 353 and the south- 
ern lobe of M 84. 

In all cases, the iso-RM contours are straight and perpendicu- 
lar to the major axes of the lobes to a very good approximation: the 
very straight and well-defined bands in the eastern lobes of both 
0206+35 (Fig.[2ji) and 3C353 (Fig.|2jl) are particularly striking. 
The entire area of M 84 appears to be covered by a banded structure, 
while in the central parts of 0206+35 and 3C 270 and the western 
lobe of 3C 353, regions of isotropic and random RM fluctuations 
are also present. 

We also derived profiles of (RM) along the radio axis of each 
source, averaging over boxes a few beamwidths long (parallel to 
the axes), but extended perpendicular to them to cover the entire 
width of the source. The boxes are all large enough to contain 
many independent points. The profiles are shown in Fig.|5] For 
each radio galaxy, we also plot an estimate of the Galactic con- 
tribution to the RM derived from a weighted mean of the integrated 
RM's fo r non-cluster radio sources wi thin a surrounding area of 
lOdeg 2 dSimard-Normandin et alj["l98ll) . In all cases, both posi- 
tive and negative fluctuations with respect to the Galactic value are 
present. 

In 0206+35 (Fig.|2ji), the largest-amplitude bands are in the 
outer parts of the lobes, with a possible low-level band just to the 
NW of the core. The most prominent band (with the most negative 
RM values) is in the eastern (receding) lobe, about 15 kpc from the 
core (Fig.|5ji). Its amplitude with respect to the Galactic value is 
about 40radm~ 2 . This band must be associated with a strong or- 
dered magnetic field component along the line of sight. If corrected 
for the Galactic contribution, the two adjacent bands in the eastern 
lobe would have RM with opposite signs and the field component 
along the line of sight must therefore reverse. 

M84 (Fig.[2})) displays an ordered RM pattern across the 
whole source, with two wide bands of opposite sign having the 
highest absolute RM values. T here is also an abrupt change of sign 
across the radio core (see also lLaing &~Br idle 1987). The negative 
band in the northern lobe (associated with the approaching jet) has 
a larger amplitude with respect to the Galactic value than the cor- 
responding (positive) feature in the southern lobe (Fig.[5j:). 

3C 270 (Fig.[2j;) shows two large bands: one on the front end 
of the eastern lobe, the other in the middle of the western lobe. 
The bands have opposite signs and contain the extreme positive 
and negative values of the observed RM. The peak positive value is 
within the eastern band at the extreme end of the lobe (Fig.|5^). 
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The RM structure of 3C 353 (Fig.|2}l) is highly asymmetric. 
The eastern lobe shows a strong pattern, made up of four bands, 
with very straight iso-RM contours which are almost exactly per- 
pendicular to the source axis. As in 0206+35, adjacent bands have 
RM with opposite signs once corrected for the Galactic contribu- 
tion (Fig.[5g). In contrast, the RM distribution in the western lobe 
shows no sign of any banded structure, and is consistent with ran- 
dom fluctuations superimposed on an almost linear profile. It seems 
very likely that the differences in RM morphology and axial ratio 
are both related to the external density gradient (Fig.[T]l). 

In Table [4] the relevant geometrical features (size, distance 
from the radio core, (RM) ) for the RM bands are listed. 



5 DEPOLARIZATION 

In this section, we use "depolarization" in its conventional sense 
to mean "decrease of degree of polarization with increasing wave- 
length " and define DP = pi.4GHz/p4.9GHz. Using the Faraday 
code dMurgia etal]|2004) . we produced images of Burn law k 
by weighted least-squares fitting to lnp(A) as a function of A 4 
(Eq.[3}. Only data with signal-to-noise ratio >4 in P at each fre- 
quency were included in the fits. The Burn law k images were pro- 
duced with the same angular resolutions as the RM images. The 
1.65arcsec resolution Bum law k maps for M84 and 3C270 are 
consistent with the low-resolution ones, but add no additional detail 
and are quite noisy. This could lead to signi ficantly biased e stimates 
for the mean values of k over large areas dLaing et alj2008[) . There- 
fore, as for the RM maps, we used only the Bum law k images at 
low resolution for these two sources. 

The Bum law k maps are shown in Fig.[3] All of the sources 
show low average values of k (i.e. slight depolarization), suggest- 
ing little RM power on small scales. With the possible exception 
of the narrow filaments of high k in the eastern lobe of 3C 353 
(which might result from partially resolved RM gradients at the 
band edges), none of the images show any obvious structure re- 
lated to the RM bands. For each source, we have also compared the 
RM and Burn law k values derived by averaging over many small 
boxes covering the emission, and we find no correlation between 
them. 

We also derived profiles of k (Fig.[5p, d, f and h) with the same 
sets of boxes as for the RM profiles in the same Figure. These con- 
firm that the values of k measured in the centres of the RM bands 
are always low, but that there is little evidence for any detailed cor- 
relation. 

The signal-to-noise ratio for 0206+35 is relatively low com- 
pared with that of the other three sources, particularly at 4.9 GHz 
(we need to use a small beam to resolve the bands), and this is 
reflected in the high proportion of blanked pixels on the k image. 
The most obvious feature of this image (Fig.[3^), an apparent differ- 
ence in mean k between the high-brightness jets (less depolarized) 
and the surrounding emission, is likely to be an artefact caused by 
our blanking strategy: points where the polarized signal is low at 
4.9 GHz are blanked preferentially, so the remainder show artifi- 
cially high polarization at this frequency. For the same reason, the 
apparent minimum in k at the centre of the deep, negative RM band 
(Fig. [5^ and b) is probably not significant. The averaged values of 
k for 0206+35 are already very low, however, and are likely to be 
slightly overestimated, so residual RM fluctuations on scales below 
the 1.2-arcsec beamwidth must be very small. 

M 84 shows one localised area of very strong depolarization 
(k ~500rad 2 m -4 , corresponding to DP = 0.38) at the base of 



the southern jet (Fig.[3j3). There is no corresponding feature in the 
RM image (Fig. [2p). The depolarization is likely to be associated 
with one of the shells of compressed gas visible in the Chandra 
image (Fig. [TJj), implying significant magnetization with inhomo- 
geneous field and/or density structure on scales much smaller than 
the beamwidth, apparently independent of the larger-scale field re- 
sponsible for the RM bands. This picture is supported by the good 
spatial coincidence of the high k region with a shell of compressed 
gas, as illustrated in the overlay of the 4.5 arcsec Burn law k image 
on the contours of the Chandra data (Fig.[6fa)). Cooler gas associ- 
ated with the emission-line disk might also be responsible, but there 
is no evidence for spatial c oincidence between enhanced depolar- 
izatio n and Ha emission ^Hansen. Norgaard-Nielsen & Jorgensenl 
1985). Despite the complex morphology of the X-ray emission 
around M84, its k profile is very symmetrical, with the highest 
values at the centre (Fig.|5jd)). 

3C 270 also shows areas of very strong depolarization 
(k -~550rad 2 m -4 , corresponding to DP = 0.35) close to the core 
and surrounding the inner and northern parts of both the radio 
lobes. As for M 84, the areas of high k are coincident with ridges 
in the X-ray emission which form the boundaries of the cavity sur- 
rounding the lobes (Fig.[6jb)). The in ner parts of this X-r ay struc- 
ture are described in more detail by Worr all et al. whose 
recent high-resolution Chandra image clearly reveals "wedges" of 
low brightness surrounding the radio jets. As in M 84, the most 
likely explanation is that a shell of denser gas immediately sur- 
rounding the radio lobes is magnetized, with significant fluctuations 
of field strength and density on scales smaller than our 5-arcsec 
beam, uncorrelated with the RM bands. The k profile of 3C 270 
(Fig-Uff)) is ver Y symmetrical, suggesting that the magnetic-field 
and density distributions are also symmetrical and consistent with 
an orientation close to the plane of the sky. The largest values of k 
are observed in the centre, coincident with the features noted ear- 
lier and with the bulk of the X-ray emission (the high k values in 
the two outermost bins have low signal-to-noise and are not signif- 
icant). 

In the Burn law k image of 3C 353, there is evidence for a 
straight and knotty region of high depolarization ~20 kpc long and 
extending westwards from the core. This region does not appear 
to be related either to the jets or to any other radio feature. As in 
M 84 and 3C 270, the RM appears quite smooth over the area show- 
ing high depolarization, again suggesting that there are two scales 
of structure, one much smaller than the beam, but producing zero 
mean RM and the other very well resolved. In 3C353, there is as 
yet no evidence for hot or cool ionized gas associated with the en- 
hanced depolarization (contamination from the very bright nuclear 
X-ray emission affects an area of 1 arcmi n radius around the core; 
llwasawa et alfeOOdlGoodger et alj2()08h . 

The k profile of 3C353 (Fig.|5jh)) shows a marked asymme- 
try, with much higher values in the East. This is in the same sense as 
the difference of RM fluctuation amplitudes (Fig.|5jg)) and is also 
consistent with the eastern lobe being embedded in higher-density 
gas. The relatively high values of k within 20 kpc of the nucleus in 
the Western lobe are due primarily to the discrete region identified 
earlier. 



6 ROTATION MEASURE STRUCTURE FUNCTIONS 

We calculated RM structure function for discrete regions of the 
sources where the RM fluctuations appear to be isotropic and ran- 
dom and for which we expect the spatial variations of foreground 
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Table 4. Properties of the RM bands: Col. 1 source name; Cols. 2&3: overall (RM) and ct rm ; Col. 4: Galactic (RM) ; Col. 5: (RM) for each band; Col. 
6: distance of the band midpoint from the radio core (positive distances are in the western direction for all sources, except for M 84, where they are in the 
northern direction); Col. 7: width of the band; Col. 8: maximum band amplitude. 
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thermal gas density, rms magnetic field strength and path length to 
be reasonably small. These are: the inner 26 arcsec of the receding 
(Eastern) lobe of 0206+35, the inner 100 arcsec of 3C270 and the 
inner 40 arcsec of the western lobe of 3C353. The selected areas 
of 0206+35 and 3C 270 are both within the core radii of the larger- 
scale beta models that describe the group-scale X-ray emission and 
the galaxy-scale components are too small to affect the RM statis- 
tics significantly (Table[2]l. In 3C 353, the selected area was chosen 
to be small compared with the scale of X-ray variations seen in 
Fig. UJd). In all three cases, the foreground fluctuations should be 
fairly homogeneous. There are no suitable regions in M 84, which 
is entirely covered by the banded RM pattern. 

The structure functions, corrected for uncorr elated random 
noise by subtracting 2cr 2 oise dSimonetti et alJ[T984l) . are shown in 
Fig.|7] All of the observed structure functions correspond to power 
spectra of approximately power-law form over all or most of the 
range of spatial frequencies we sample. We initially assumed that 
the power spectrum was described by Eq[5]with no high-frequency 
cut-off (/ m ax — > oo) and made least-squares fits to the structure 
functions, weighted by errors derived from multiple realizations of 
the power spectrum on the observing grid, a s described in detail by 
lLaing et alj <2008b and lGuidetti et all koid) . 

The best-fitting slopes q and amplitudes Co are given in Ta- 
ble[5] All of the fitted power spectra are quite flat and have low 
amplitudes, implying that there is little power in the isotropic and 
random component of rotation measure. Indeed, the amplitudes of 
the largest-scale RM fluctuations sampled in this analysis is a few 
times less than that of the bands (see Tables|4]and[5j. This suggests 
that the field responsible for the bands is stronger as well as more 
ordered than that responsible for the isotropic fluctuations. 

The structure functions for 0206+35 and 3C 353 rise mono- 
tonically, indicating that the outer scale for the random fluctua- 
tions must be larger than the maximum separations we sample. 
For 3C270, the structure function levels out at r± « 100 arcsec 
(15 kpc; Fig. [7}l). This could be the outer scale of the field fluc- 
tuations, but a better understanding of the geometry and external 
density distribution would be needed before we could rule out the 



effects of large-scal e variations in path length or field strength (cf. 
iGuidetti et alfeoiol) . 

In order to constrain RM structure on spatial scales below the 
beamwidth, we estimated the depolarization as described in Sec- 
tion^ The fitted k values are listed in Tablef5] We stress that these 
values refer only to areas with isotropic fluctuations, and cannot 
usefully be compared with the integrated depolarizations quoted in 
in Fig.|2] 

For M 84, using the Burn law k analysis and assuming that 
variation of Faraday rotation across the 1.65-arcsec beam causes 
the residual depolarization, we find that A m i n < 0.1 kpc for any rea- 
sonable RM power spectrum. 



7 ROTATION-MEASURE BANDS FROM COMPRESSION 

It is clear from the fact that the observed RM bands are perpen- 
dicular to the lobe axes that they must be associated with an in- 
teraction between an expanding radio source and the gas imme- 
diately surrounding it. One inevitable mechanism is enhancement 
of field and density by the shock or compression wave surrounding 
the source0The implication of the presence of cavities in the X-ray 
gas distribution coincident with the radio lobes is that the sources 
are interacting strongl y with the thermal gas, dis placing rather than 
mixing with it (see iMcNamara & Nulsenll2007l for a review). For 
the sources in the presen t paper, the X-ray o bservations of M 84 
dFinoguenov et al. ll2008l FigQJi) and 3C 270 dCroston et ai]|2008l 
FigOl) show cavities and arcs of enhanced brightness, correspond- 
ing to shells of compressed gas bounded by weak shocks. The 
strength of any pre-existing field in the IGM, which will be frozen 
into the gas, will also be enhanced in the shells. We therefore ex- 
pect a significant enhancement in RM. A more extreme example of 
this effect will occur if the expansion of the radio source is highly 
supersonic, in which case there will be a strong bow-shock ahead of 

2 An alternative mechanism is the g eneration of non-linear surface waves 
(Bicknell. Cameron & Gi ngold 1990). It is unlikely that this can produce 
large-scale bands, for the reasons given in SectionRHl 
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Table 5. Power spectrum parameters for the individual sub-regions. Col. 1: source name; Col. 2: angular resolution; Col. 3: slope q: Col. 4: amplitude Co; 
Col. 5: minimum scale; Col. 6: amplitude of the large scale isotropic component; Col. 7: observed mean k; Col. 8: predicted mean k. The power spectrum has 
not been computed for M 84 (see Section|6j. 
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Figure 4. Plots of E- vector position angle against A 2 at representative points of the 1.2 arcsec RM map of 0206+35. Fits to the relation VP (A) = ^o+RMA 2 
are shown. The values of RM are given in the individual panels. 



the lobe, behind which both the density and the field become much 
higher. Regardless of the strength of the shock, the field is modified 
so that only the component in the plane of the shock is amplified 
and the post-shock field tends become ordered parallel to the shock 
surface. 

The evidence so far suggests that shocks around radio sources 
of both FR classe s are generally wea k (e.g. iForman et al. 1120051 
Wilso n et al. I l200d iNulsen et al.ll2005h . There are only two exam- 
ples in which highly supersonic expansion has been inferre d: the 
southern lobe of C entaurusA (M ~ 8; iKraft et al.l 120031) and 
NGC 3801 (M w 4; ICroston et alj2007l) . There is no evidence that 
the sources described in the present paper are significantly over- 
pressured compared with the surrounding IGM (indeed, the syn- 
chrotron minimum pressure is systematically lower than the ther- 
mal pressure of the IGM; Table |2j. The sideways expansion of the 
lobes is therefore unlikely to be highly supersonic. The shock Mach 
number estimated for all the sources from ram pressure balance in 
the forward direction is also ~ 1.3. This estimat e is consistent with 
that for M 84 made by iFinoguenov et aT I d2006h and also with the 



lack of detection of strong shocks in the X-ray data for the other 
sources. 

In this section, we investigate how the RM could be affected 
by compression. We consider a deliberately oversimplified picture 
in which the radio source expands into an IGM with an initially 
uniform magnetic field, B. This is the most favourable situation 
for the generation of large-scale, anisotropic RM structures: in re- 
ality, the pre-existing field is likely to be highly disordered, or even 
isotropic, because of turbulence in the thermal gas. We stress that 
we have not tried to generate a self-consistent model for the mag- 
netic field and thermal density, but rather to illustrate the generic 
effects of compression on the RM structure. 

In this model the radio lobe is an ellipsoid with its major axis 
along the jet and is surrounded by a spherical shell of compressed 
material. This shell is centred at the mid-point of the lobe (Fig.[8j( 
and has a stand-off distance equal to 1/3 of the lobe semi-major 
axis at the leading edge (the radius of the spherical compression 
is therefore equal to 4/3 of the lobe semi-major axis). In the com- 
pressed region, the thermal density and the magnetic field compo- 
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Figure 5. Profiles of mean and rms RM and Burn law k along the radio axis. The profiles have been derived by averaging over boxes perpendicular to the radio 
axis of length 5 kpc for all sources except M 84, for which the box length is 2 kpc. The orange (dashed) and violet (long-dashed) lines represent the Galactic 
RM contribution and the (RM) of the lobes, respectively. All unblanked pixels are included. 
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Figure 7. Plots of the RM structure functions for the isotropic sub-regions of 0206+35, 3C 353 and 3C 270 described in the text. The horizontal bars represent 
the bin width and the crosses the centroids for data included in the bins. The red lines are the predictions for power law power spectra, including the effects of 
the convolving beam. The vertical error bars are the rms variations for the structure functions derived for multiple realizations of the data. The fitted structure 
functions for 3C 270 are derived for the same power spectrum parameters. 
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Figure 8. Amplification and sweeping up of the magnetic field lines inside 
the compression region defined by the projected circular green and dashed 
line. The driving expansion is along the Z-axis and the z-axis represents a 
generic line-of- sight. 



nent in the plane of the spherical compression are amplified by the 
same factor, because of flux-freezing. We use a coordinate system 
xyz centred at the lobe mid-point, with the 2-axis along the line 
of sight, so x and y are in the plane of the sky. The radial unit 
vector is r=(x,y, z). ~B=(B X , B y , B z ) and ~B'=(B' X , B' y , B' z ) are 
respectively the pre- and post-shock magnetic-field vectors. Then, 
we consider a coordinate system XYZ still centred at the lobe mid- 
point, but rotated with respect to the xyz system by the angle 6 
about the y (Y) axis, so that Z is aligned with the major axis of 
the lobe. With this choice, 6 is the inclination of the source with 
respect to the line-of-sight (Fig.[8j. 

After a spherical compression, the thermal density and field 
satisfy the equations: 

n e = jn e 
B' ± = B ± = (B r)r 

B|, = 7 B|, = 7 [B-(B-r)f] (5) 

where n c , Bi and Bii represent the initial thermal density and the 
components of the field perpendicular and parallel to the compres- 
sion surface. The same symbols with primes stand for post-shock 
quantities and 7 is the compression factor. The total compressed 
magnetic field is: 



B' = B|, + B' ± = 7B + (1 - 7 )(B ■ r)f 



(6) 



The field strength after compression depends on the the angle 
between the compression surface and the initial field. Maximum 
amplification occurs for a field which is parallel to the surface, 
whereas a perpendicular field remains unchanged. The post-shock 
field component along the line-of-sight becomes: 



B' z = B' ■ z = 7B ■ z + (1 - 7)(B ■ r)(r ■ 



(7) 



We assumed that the compression factor, 7 = j(Z), is a func- 
tion of distance Z along the source axis from the centre of the radio 
lobe, decreasing monotonically from a maximum value 7 max at the 
leading edge to a constant value from the centre of the lobe as far as 
the core. We investigated values of 7 max in the range 1 .5 - 4 (7 = 4 
corresponds to the asymptotic value for a strong shock). Given that 



there is no evidence for strong shocks in the X-ray data for any of 
our sources, we have typically assumed that the compression factor 
is 7max = 3 at the front end of the lobe, decreasing to 1.2 at the 
lobe mid-point and thereafter remaining constant as far as the core. 
A maximum compression factor of 3 is consistent with the tran- 
sonic Mach numbers M — 1.3 estimated from ram-pressure bal- 
ance for all of the sources and this choice is also motivated by the 
X-ray data of M 84, from which there is evidence for a compression 
ratio « 3 between the shells and their surroundings (Section |2!4l . 

We produced synthetic RM images for different combinations 
of source inclination and direction of the pre-existing uniform field, 
by integrating the expression 



RM syn = 



n'B'dz. 



(8) 



2.4 x 



numerically. We assumed a constant value of n e 
10 -3 cm -3 for the density of the pre-shock material (the central 
value for the group gas associated with 0206+35; Table [2}, a lobe 
semi-major axis of 21 kpc (also appropriate for 0206+35) and an 
initial field strength of 1 /iG. The integration limits were defined by 
the surface of the radio lobe and the compression surface. This is 
equivalent to assuming that there is no thermal gas within the ra- 
dio lobe, consistent with the picture suggested by our inference of 
foreground Faraday rotation and the existence of X-ray cavities and 
that Faraday rotation from uncompressed gas is negligible. 

As an example, Fig.[9]shows the effects of compression on the 
RM for the receding lobe of a source inclined by 40° to the line 
of sight. The initial field is pointing towards us with an inclination 
of 60° with respect to the line-of-sight; its projection on the plane 
of the sky makes an angle of 30° with the a;-axis. Fig.[9ja) displays 
the RM produced without compression (7(Z) = 1 everywhere): the 
RM structures are due only to differences in path length across the 
lobe. Fig.|9lb) shows the consequence of adding a modest compres- 
sion of 7 m ax = 1.5: structures similar to bands are generated at the 
front end of the lobe and the range of the RM values is increased. 
Fig.|9jc) illustrates the RM produced in case of the strongest pos- 
sible compression, 7 ma x = 4: the RM structure is essentially the 
same as in Fig.[9jb), with a much larger range. 

This very simple example shows that RM bands with ampli- 
tudes consistent with those observed can plausibly be produced 
even by weak shocks in the IGM, but the iso-RM contours are 
neither straight, nor orthogonal to the lobe axis and there are no 
reversals. These constraints require specific initial conditions, as il- 
lustrated in Fig.[l0] where we show the RM for a lobe in the plane 
of the sky. We considered three initial field configurations: along 
the line-of-sight (Fig.UOb), in the plane of the sky and parallel to 
the lobe axis (Fig.llOb) and in the plane of the sky, but inclined by 
45° to the lobe axis (Fig.UOfc). The case closest to reproducing the 
observations is that displayed in Fig.llOfb), in which reversals and 
well defined and straight bands perpendicular to the jet axis are pro- 
duced for both of the lobes. In Fig.UOfa). the structures are curved, 
while in Fig JlOf c) the bands are perpendicular to the initial field 
direction, and therefore inclined with respect to the lobe axis. 

For a source inclined by 40° to the line of sight, we found 
structures similar to the observed bands only with an initial field in 
the plane of the sky and parallel to the axis in projection (Figsll lb 
and b; note that the synthetic RM images in this example have been 
made for each lobe separately, neglecting superposition). 

We can summarize the results of the spherical pure compres- 
sion model as follows. 
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Figure 9. Panel (a): synthetic RM for a receding lobe inclined by 40° to the line-of-sight and embedded in a uniform field inclined by 60° to the line-of-sight 
and by and 30° to the lobe axis on the plane of the sky. Panels (b) and (c): as (a) but with a weak compression (7 ma x = 1-5) and a strong compression 
(7max = 4) of the density and field. The crosses and the arrows indicate the position of the core and the lobe advance direction, respectively. 



(i) An initial field with a component along the line-of-sight does 
not generate straight bands. 

(ii) The bands are orthogonal to the direction of the initial field 
projected on the plane of the sky, so bands perpendicular to the lobe 
axis are only obtained with an initial field aligned with the radio jet 
in projection. 

(iii) The path length (determined by the precise shape of the ra- 
dio lobes) has a second-order effect on the RM distribution (com- 
pare Figs.|9ji and b). 

Thus, a simple compression model can generate bands with 
amplitudes similar to those observed but reproducing their geome- 
try requires implausibly special initial conditions, as we discuss in 
the next section. 



8 ROTATION-MEASURE BANDS FROM A DRAPED 
MAGNETIC FIELD 

8.1 General considerations 

That the pre-existing field is uniform, close to the plane of the sky 
and aligned with the source axis in projection is implausible for 
obvious reasons: 

(i) the pre-existing field cannot know anything about either the 
radio-source geometry or our line-of-sight and 

(ii) observations of Faraday rotation in other sources and the 
theoretical inference of turbulence in the IGM both require disor- 
dered initial fields. 

This suggests that the magnetic field must be aligned by rather 
than with the expansion of the radio source. Indeed, the field con- 
figurations which generate straigh t bands look qualitatively like 
the "draping" model proposed by iDursi & Pfrommed ( 1200 81) . for 
some angles to the line-of-sight. The analysis of Sections[2][5]and 
[6] suggests that the magnetic fields causing the RM bands are 
well-ordered, consistent with a stretching of the initial field that 
has erased much of the small-scale structure while amplifying the 
large-scale component. We next attempt to constrain the geometry 
of the resulting "draped" field. 

8.2 Axisymmetric draped magnetic fields 

A proper c alculation of the RM from a draped magnetic-field con- 
figuration dDursi & Pfromm er 2008) is outside the scope of this 
paper, but we can start to understand the field geometry using 



some simple approximations in which the field lines are stretched 
along the source axis. We assumed initially that the field around 
the lobe is axisymmetric, with components along and radially out- 
wards from the source axis, so that the RM pattern is independent of 
rotation about the axis. It is important to stress that such an axisym- 
metric field is not physical, as it requires a monopole and unnatu- 
ral reversals, it is nevertheless a useful benchmark for features of 
the field geometry that are needed to account for the observed RM 
structure. We first considered field lines which are parabolae with 
a common vertex on the axis ahead of the lobe. For field strength 
and density both decreasing away from the vertex, we found that 
RM structures, with iso-contours similar to arcs, rather than bands, 
were generated only for the approaching lobe of an inclined source. 
Such anisotropic RM structures were not produced in the reced- 
ing lobes, nor for sources in the plane of the sky. Indeed, in order 
to generate any narrow, transverse RM structures such as arcs or 
bands, the line-of-sight must pass through a foreground region in 
which the field lines show significant curvature, which occurs only 
for an approaching lobe in case of a parabolic field geometry. This 
suggested that we should consider field lines which are families of 
ellipses centred on the lobe, again with field strength and density 
decreasing away from the leading edge. This indeed produced RM 
structures in both lobes for any inclination, but the iso-RM contours 
were arcs, not straight lines. Because of the non-physical nature of 
these axisymmetric field models, the resulting RM images are de- 
liberately not shown in this paper. In order to quantify the depar- 
tures from straightness of the iso-RM contours, we measured the 
ratio of the predicted RM values at the centre and edge of the lobe 
at constant Y, at different distances along the source axis, X, for 
both of the example axisymmetric field models. The ratio, which 
is 1 for perfectly straight bands, varies from 2 to 3 in both cases, 
depending on distance from the nucleus. This happens because the 
variations in line-of-sight field strength and density do not com- 
pensate accurately for changes in path length. We believe that this 
problem is generic to any axisymmetric field configuration. 
The results of this section suggest that the field configuration re- 
quired to generate straight RM bands perpendicular to the projected 
lobe advance direction has systematic curvature in the field lines (in 
order to produce a modulation in RM) without a significant depen- 
dence on azimuthal angle around the source axis). We therefore 
investigated a structure in which elliptical field lines are wrapped 
around the front of the lobe, but in a two-dimensional rather than a 
three-dimensional configuration. 
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Figure 10. Synthetic RM images for a lobe in the plane of the sky with an ambient field which is uniform before compression. In panels (a) and (b), the field is 
along the axis of the lobe in projection and inclined to the line-of-sight by (a) 45° and (b) 90° . In panel (c), the field is in the plane of the sky and misaligned 
by 45° with respect to the lobe axis. The crosses and the arrows represent the radio core position and the lobe advance direction, respectively. 





Figure 11. Synthetic RM images of the receding (a) and approaching lobe (b) of a source inclined by 40° with respect to the line-of-sight, in the case of a 
uniform ambient field. The field is aligned with the radio jets in projection and inclined with respect to the line-of-sight by 90° . The crosses and the arrows 
represent the radio core position and the lobe advance direction, respectively. 



8.3 A two-dimensional draped magnetic field 

We considered a field with a two-dimensional geometry, in which 
the field lines are families of ellipses in planes of constant Y, as 
sketched in Fig.[l2] The field structure is then independent of Y. 
The limits of integration are given by the lobe surface and an ellipse 
whose major axis is 4/3 of that of the lobe. Three example RM im- 
ages are shown in Fig.[l3] The assumed density (2.4 x 10 -3 cm -3 ) 
and magnetic field (1 /iG) were the same as the pre-shock values 
for the compression model of Section [7] and the lobe semi-major 
axis was again 21 kpc. Since the effect of path length is very small 
(Section[7), the RM is also independent of Y to a good approxima- 
tion. The combination of elliptical field lines and invariance with Y 
allows us to produce straight RM bands perpendicular to the pro- 
jected lobe axis for any source inclination. Furthermore, this model 
generates more significant reversals of the RM (e.g. Fig.ll3b) than 
those obtainable with pure compression (e.g. Fig.ll lb). 

We conclude that a field model of this generic type represents 
the simplest way to produce RM bands with the observed character- 
istics in a way that does not require improbable initial conditions. 
The invariance of the field with the Y coordinate is an essential 
point of this model, suggesting that the physical process respon- 
sible for the draping and stretching of the field lines must act on 
scales larger than the radio lobes in the Y direction. 

8.4 RM reversals 

The two-dimensional draped field illustrated in the previous section 
reproduces the geometry of the observed RM bands very well, but 
can only generate a single reversal, which must be very close to the 



front end of the approaching lobe, where the elliptical field lines 
bend most rapidly. We observe a prominent reversal in the receding 
lobe of 0206+35 (Fig. [2^) and multiple reversals across the eastern 
lobe of 3C353 and in M84 (Fig. [2jj and d). The simplest way to 
reproduce these is to assume that the draped field also has rever- 
sals, presumably originating from a more complex initial field in 
the IGM. One realization of such a field configuration would be in 
the form of multiple toroidal eddies with radii smaller than the lobe 
size, as sketched in Fig.[l4] Whatever the precise field geometry, 
we stress that the straightness of the observed multiple bands again 
requires a two-dimensional structure, with little dependence on Y. 



9 DISCUSSION 

9.1 Where do bands occur? 

The majority of published RM images of radio galaxies do not 
show bands or other kind of anisotropic structure, but are charac- 
terized by isotropic and ra ndom RM distributions (e.g. lLaing et al.l 
2008; Guid etti et alj|2010f) . On the other hand, we have presented 
observations of RM bands in four radio galaxies embedded in dif- 
ferent environments and with a range of jet inclinations with respect 
to the line-of-sight. Our sources are not drawn from a complete 
sample, so any quantitative estimate of the incidence of bands is 
premature, but we can draw some preliminary conclusions. 

The simple two-dimensional draped-field model developed in 
Section [831 only generates RM bands when the line-of-sight inter- 
cepts the volume containing elliptical field lines, which happens for 
a restricted range of rotation around the source axis. At other ori- 
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Figure 12. Geometry of the two-dimensional draped model for the magnetic field, described in Sec. 18.31 Panel (a): seen in the plane of the direction of the 
ambient magnetic field. Panel (b): seen in the plane perpendicular to the ambient field; only the innermost family of ellipses is plotted. For clarity, field lines 
behind the lobe are omitted. The colour of the curves indicates the strength of the field, decreasing from red to yellow. The crosses represent the radio core 
position. The coordinate system is the same as in Fig.[8]and the line of sight is in the X — Z plane. 
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Figure 13. Synthetic RM images produced for a two-dimensional draped model as described in Section 18.31 The field lines are a family of ellipses draped 
around the lobe as in Fig. 1121 Panel (a): lobe inclined at 90° to the line-of-sight. (b) receding and (c) approaching lobes of a source inclined by 40° to the 
line-of-sight. The crosses and the arrows represent the radio core position and the lobe advance direction. 



entations, the RM from this field configuration will be small and 
the observed RM may well be dominated by material at larger dis- 
tances which has not been affected by the radio source. We there- 
fore expect a minority of sources with this type of field structure 
to show RM bands and the remainder to have weaker, and prob- 
ably isotropic, RM fluctuations. In contrast, the thr ee-dimensional 
draped field model proposed bv lDursi & Pfrommer| J2008h predicts 
RM bands parallel to the source axis for a significant range of view- 
ing directions: these have not (yet) been observed. 



9.2 RM bands in other sources 

The fact that RM bands have so far been observed only in a few 
radio sources may be a selection effect: much RM analysis has been 
carried out for galaxy clusters , in which most of the sources are 
tailed (e.g. Blanton et al. 2003). With a few exceptions like Cyg A 
(see below, Section |9. 2. U , lobed FRI and FRII sources have not 
been studied in detail. 



The prominent RM bands described in the present paper oc- 
cur only in lobed radio galaxies. In contrast, well-observed ra- 
dio sources with tails and plumes seem to b e free of bands or 
anisotropic RM str ucture (e.g. 3C31, 3C 449: ITaing et alj|2008l ; 
iGuidetti et aHl2010h . Furthermore, the lobes which show bands are 
all quite round and show evidence for interaction with the surround- 
ing IGM. It is particularly striking that the bands in 3C 353 occur 
only in its eastern, rounded, lobe. The implication is that RM bands 
occur when a lobe is being actively driven by a radio jet into a re- 
gion of high IGM density. Plumes and tails, on the other hand, are 
likely to be rising buoyantly in the group or cluster and we do not 
expect significant compression, at least at large distances from the 
nucleus. 



9.2.7 CygnusA 

Cygnus A is a source in which we might expect to observed RM 
bands, by analogy with the sources discussed in the present paper: 
it has wide and round lobes and Chandra X -ray data have sho wn the 
presence of shock-heated gas and cavities JWilson et al.l2006h . RM 
bands, roug hly perpendicular to the source axis, are ind eed seen in 
both lobes jDreher et al]|l987l : ICarilli & Tavlorl 120021) , but inter- 
pretation is complicated by the larger random RM fluctuations and 
the strong depolarization in the eastern lobe. A semi-circular RM 
feature around one of the hot-spots in the western lobe has been at- 
tribut ed to compression by the bow-shock (Carilli, Perlev & Dreherl 
Il988l) . 
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Figure 14. Geometry of a two-dimensional toroidal magnetic field as de- 
scribed in Section liOl seen in the plane normal to its axis. The cross repre- 
sents the radio core position. 



9.2.2 Hydra A 

ICariili&Tavio^ ( l2002[) have claimed evidence for RM bands in the 
north ern lobe of Hydra A. The Chandra image (McNamara et al. 
2000) shows a clear cavity with sharp edges coincident with the ra- 
dio lobes and an absence of shock-heated gas, just as in our sources. 
Despite the classification as a tailed source, it may well be that there 
is significant compression of the IGM. Note, however, that the RM 
image is not well sampled close to the nucleus. 



9.2.3 3C465 

The RM ima ge of the tailed source 3C465 published by 

Eilek & OwenN2002h shows some evidence for bands, but the 
colour scale was deliberately chosen to highlight the difference 
between positive and negative values, thus making it difficult to 
see the large gradients in RM expected at band edges. The original 
RM image (Eilek, private communication) suggests that the band in 
the western tail of 3C 465 is similar to those we have identified. It 
is plausible that magnetic-field draping happens in wide-angle tail 
sources like 3C 465 as a result of bulk motion of the IGM within 
the cluster potential well, as required to bend the tails. It will be in- 
teresting to search for RM bands in other sources of this type and to 
find out whether there is any relation between the iso-RM contours 
and the flow direction of the IGM. 



9.3 Foreground isotropic field fluctuations 

The coexistence in the RM images of our sources of anisotropic 
patterns with areas of isotropic fluctuations suggests that the 
Faraday-rotating medium has at least two components: one local 
to the source, where its motion significantly affects the surround- 
ing medium, draping the field, and the other from material on group 
or cluster scales which has not felt the effects of the source. This 
raises the possibility that turbulence in the foreground Faraday ro- 
tating medium might "wash out" RM bands, thereby making them 
impossible to detect. The isotropic RM fluctuations observed across 
our sources are all described by quite flat power spectra with low 
amplitude (Table|5]l. The random small-scale structure of the field 
along the line-of-sight essentially averages out, and there is very 
little power on scales comparable with the bands. If, on the other 
hand, the isotropic field had a steeper power spectrum with sig- 



nificant power on scales similar to the bands, then its contribution 
might become dominant. 

We first produced synthetic RM images for 0206+35 including 
a random component derived from our best-fitting power spectrum 
(Table[5} in order to check that the bands remained visible. We as- 
sumed a minimum scale of 2 kpc from the depolarization analysis 
for 0206+35 (Sec.O, and a maximum scale of A max = 40 kpc, 
consistent with the continuing rise of the RM structure function at 
the largest sampled separations, which requires A max > 30 arcsec 
(~20kpc). The final synthetic RM is given by: 

RM syn = RM drap + RM icm = / n' c B' z dz + [ n c B z dz (9) 



where RM drap and RM icm are the RM due to the draped and 
isotropic fields, respectively. The terms n' c and B'^ are respec- 
tively the density and field component along the line-of-sight in 
the draped region. The integration limits of the term RM drap were 
defined by the surface of the lobe and the draped region, while that 
of the term RM lcm starts at the surface of the draped region and ex- 
tends to 3 times the core radius of the X-ray gas (Tabled . For the 
electron gas density n c outside the draped region we assumed the 
beta-model profile of 0206+35 (Table[2} and for th e field strength 
a radial variation of the form dGuidetti et alJl201oL and references 
therein). 



{B 2 (r)) 1/2 = B 



n e (r) 



no 



(10) 



where Bo is the rms magnetic field strength at the group centre. We 
took a draped magnetic field strength of 1.8 /iG, in order to match 
the amplitudes for the RM bands in both lobes of 0206+35, and 
assumed the same value for Bq . 

Example RM images, shown in Figs.|15[a) and (b), should 
be compared with those for the draped field alone (Figs.ll3b and 
c, scaled up by a factor of 1.8 to account for the difference in 
field strength) and with the observations (Fig(2^) after correction 
for Galactic foreground (Table|4j. The model is self-consistent: the 
flat power spectrum found for 0206+35 does not give coherent RM 
structure which could interfere with the RM bands, which are still 
visible. 

We then replaced the isotropic component with one having 
a Kolmogorov power spectrum (q — 11/3). We assumed identi- 
cal minimum and maximum scales (2 and 40 kpc) as for the pre- 
vious power spectrum and took the same central field strength 
(£>o=l-8/iG) and radial variation (Eq.llOt. Example realizations are 
shown in Fig. |15f c) and (d). The bands are essentially invisible in 
the presence of foreground RM fluctuations with a steep power 
spectrum out to scales larger than their widths. 

It may therefore be that the RM bands in our sources are 
especially prominent because the power spectra for the isotropic 
RM fluctuations have unusually low amplitudes and flat slopes. We 
have established these parameters directly for 0206+35, 3C 270 and 
3C 353; M 84 is only 14 kpc in size and is located far from the core 
of the Virgo cluster, so it is plausible that the cluster contribution to 
its RM is small and constant. 

We conclude that the detection of RM bands could be influ- 
enced by the relative amplitude and scale of the fluctuations of the 
isotropic and random RM component compared with that from any 
draped field, and that significant numbers of banded RM structures 
could be masked by isotropic components with steep power spec- 
tra. 
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Figure 15. Synthetic RM images of the lobes of 0206+35, produced by the sum of a draped field local to the source and an isotropic, random field spread 
through the group volume with the best-fitting power spectrum found in Sec.[6] (panels (a) and (b)) and with a Kolmogorov power spectrum (panels (c) and 
(d)). Both power spectra have the same high and low-frequency cut offs and the central magnetic field strengths are also identical. The crosses and the arrows 
represent the radio core position and the lobe advance direction. 



9.4 Asymmetries in RM bands 

The well-established correlation between RM variance and/or de- 
polarization and jet sidedness observed in FRI and FRII radio 
sources is interpreted as an orientation effect: the lobe contain- 
ing the brighter jet is on t he near side, and is seen through less 
magnetoionic material (e.g. lLain 3 19881 : iMorgantietaO 19971) . For 
sources showing RM bands, it is interesting to ask whether the 
asymmetry is due to the bands or just to the isotropic component. 

In 0206+35, whose jets are inclined by ~ 40° to the line-of- 
sight, the large negative band on the receding side has the highest 
RM (Fig.[5](. This might suggest that the RM asymmetry is due to 
the bands, and therefore to the local draped field. Unfortunately, 
0206+35 is the only source displaying this kind of asymmetry. The 
other "inclined" source, M 84, (9 ~ 60° ) does not show such asym- 
metry: on the contrary the RM amplitude is quite symmetrical. In 
this case, however, the relation between the (well-constrained) in- 
clination of the inner jets, and that of the lobes could be compli- 
cated: both j ets bend by m 90° at distances of about 50arcsec from 
the nucleus jLaing et alj|201 l[), so that we cannot establish the real 
orientation of the lobes with respect to the plane of the sky. The 
low values of the jet/counter-jet ratios in 3C 270 and 3C 353 sug- 
gest that their axes are close to the plane of the sky, so that little 
orientation-dependent RM asymmetry would be expected. Indeed, 
the lobes of 3C 270 show similar RM amplitudes, while the large 
asymmetry of RM profile of 3C 353 is almost certainly due to a 
higher column density of thermal gas in front of the eastern lobe. 
Within our small sample, there is therefore no convincing evidence 
for higher RM amplitudes in the bands on the receding side, but 
neither can such an effect be ruled out. 

In models in which an ordered field is draped around the radio 



lobes, the magnitude of any RM asymm etry depends on t he field 
geometry as well as the path length (cf. lLaing et al] |2008 for the 
isotropic case). For instance, the case illustrated in Figs|13fb) and 
(c) shows very little asymmetry even for 6= 40 ° . The presence of 
a systematic asymmetry in the banded RM component could there- 
fore be used to constrain the geometry. 



9.5 Enhanced depolarization and mixing layers 

A different m echanism for the generation of RM fl uctuations was 
suggested by iBicknell. Cameron & Gingold (1993). They argued 
that large-scale nonlinear surface waves could form on the surface 
of a radio lobe through the merging of smaller waves generated by 
Kelvin-Helmholtz instabilities and showed that RM's of roughly 
the observed magnitude would be produced if a uniform field in- 
side the lobe was advected into the mixing layer. This mechanism 
is unlikely to be able to generate large-scale bands, however: the 
predicted iso-RM contours are only straight over parts of the lobe 
which are locally flat, even in the unlikely eventuality that a coher- 
ent surface wave extends around the entire lobe. 

The idea that a mixing layer generates high Faraday rotation 
may instead be relevant to the anomalously high depolarizations 
associated with regions of compressed gas around the inner ra- 
dio lobes of M 84 and 3C 270 (Section |5] and Fig.[6]l. We have ar- 
gued that the fields responsible for the depolarization are tangled on 
small scales, since they produce depolarization without any obvi- 
ous effects on the large-scale Faraday rotation pattern. It is unclear 
whether the level of turbulence within the shells of compressed gas 
is sufficient to amplify and tangle a pre-existing field in the IGM to 
the level that it can produce the observed depolarization; a plausi- 
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ble alternative is that the field originates within the radio lobe and 
mixes with the surrounding thermal gas. 



10 CONCLUSIONS AND OUTSTANDING QUESTIONS 

In this work we have analysed and interpreted the Faraday rota- 
tion across the lobed radio galaxies 0206+35, 3C270, 3C353 and 
M 84, located in environments ranging from a poor group to one of 
the richest clusters of galaxies (the Virgo cluster). The RM images 
have been produced at resolutions ranging from 1.2 to 5.5 arcsec 
FWHM using Very Large Array data at multiple frequencies. All 
of the RM images show peculiar banded patterns across the radio 
lobes, implying that the magnetic fields responsible for the Fara- 
day rotation are anisotropic. The RM bands coexist and contrast 
with areas of patchy and random fluctuations, whose power spectra 
have been estimated using a structure-function technique. We have 
also analysed the variation of degree of polarization with wave- 
length and compared this with the predictions for the best-fitting 
RM power spectra in order to constrain the minimum scale of mag- 
netic turbulence. We have investigated the origin of the bands by 
making synthetic RM images using simple models of the interac- 
tion between radio galaxies and the surrounding medium and have 
estimated the geometry and strength of the magnetic field. 
Our results can be summarized as follows. 

(i) The lack of deviation from A 2 rotation over a wide range 
of polarization position angle and the lack of associated depolar- 
ization together suggest that a foreground Faraday screen with no 
mixing of radio-emitting and thermal electrons is responsible for 
the observed RM in the bands and elsewhere (Section[4]l. 

(ii) The dependence of the degree of polarization on wavelength 
is well fitted by a Burn law, which is also consistent with (mostly 
resolved) pure foreground rotation (Section|5]l- 

(iii) The RM bands are typically 3 - lOkpc wide and have am- 
plitudes of 10 - 50radm~ 2 (Tableg). The maximum deviations of 
RM from the Galactic values are observed at the position of the 
bands. Iso-RM contours are orthogonal to the axes of the lobes. In 
several cases, neighbouring bands have opposite signs compared 
with the Galactic value and the line-of-sight field component must 
therefore reverse between them. 

(iv) An analysis of the profiles of (RM) and depolarization 
along the source axes suggests that there is very little small-scale 
RM structure within the bands. 

(v) The lobes against which bands are seen have unusually small 
axial ratios (i.e. they appear round in projection; Fig. |2j. In one 
source (3C 353) the two lobes differ significantly in axial ratio, and 
only the rounder one shows RM bands. This lobe is on the side of 
the source for which the external gas density is higher. 

(vi) Structure function and depolarization analyses show that 
flat power-law power spectra with low amplitudes and high- 
frequency cut-offs are characteristic of the areas which show 
isotropic and random RM fluctuations, but no bands (Section|6]l. 

(vii) There is evidence for source-environment interactions, 
such as large-scale asymmetry (3C 353) cavities and shells of 
swept-up and compressed material (M 84, 3C 270) in all three 
sources for which high-resolution X-ray imaging is available. 

(viii) Areas of strong depolarization are found around the edges 
of the radio lobes close to the nuclei of 3C 270 and M 84. These are 
probably associated with shells of compressed hot gas. The absence 
of large-scale changes in Faraday rotation in these features suggests 
that the field must be tangled on small scales (Section|5j- 



(ix) The comparison of the amplitude of (RM) with that of the 
structure functions at the largest sampled separations is consistent 
with an amplification of the large scale magnetic field component 
at the position of the bands. 

(x) We produced synthetic RM images from radio lobes expand- 
ing into an ambient medium containing thermal material and mag- 
netic field, first considering a pure compression of both thermal 
density and field, and then including three- and two-dimensional 
stretching ("draping") of the field lines along the direction of the 
radio jets (Sects. |7]and[8}. Both of the mechanisms are able to gen- 
erate anisotropic RM structure. 

(xi) To reproduce the straightness of the iso-RM contours, a 
two-dimensional field structure is needed. In particular, a two- 
dimensional draped field, whose lines are geometrically described 
by a family of ellipses, and associated with compression, repro- 
duces the RM bands routinely for any inclination of the sources to 
the line-of-sight (Sec. 18. 31 ). Moreover, it might explain the high RM 
amplitude and low depolarization observed within the bands. 

(xii) The invariance of the magnetic field along the axis per- 
pendicular to the forward expansion of the lobe suggests that the 
physical process responsible for the draping and stretching of the 
magnetic field must act on scales larger than the lobe itself in this 
direction. We cannot yet constrain the scale size along the line of 
sight. 

(xiii) In order to create RM bands with multiple reversals, 
more complex field geometries such as two-dimensional eddies are 
needed (Section[8]4](. 

(xiv) We have interpreted the observed RM's as due to two mag- 
netic field components: one draped around the radio lobes to pro- 
duce the RM bands, the other turbulent, spread throughout the sur- 
rounding medium, unaffected by the radio source and responsible 
for the isotropic and random RM fluctuations (Section 19.31 ). We 
tested this model for 0206+35, assuming a typical variation of field 
strength with radius in the group atmosphere, and found that a mag- 
netic field with central strength of Bq = 1.8 /iG reproduced the 
RM range quite well in both lobes. 

(xv) We have suggested two reasons for the low rate of detection 
of bands in published RM images: our line of sight will only inter- 
cept a draped field structure in a minority of cases and rotation by 
a foreground turbulent field with significant power on large scales 
may mask any banded RM structure. 

Our results therefore suggest a more complex picture of the 
magnetoionic environments of radio galaxies than was apparent 
from earlier work. We find three distinct types of magnetic-field 
structure: an isotropic component with large-scale fluctuations, 
plausibly associated with the undisturbed intergalactic medium; a 
well-ordered field draped around the leading edges of the radio 
lobes and a field with small-scale fluctuations in the shells of com- 
pressed gas surrounding the inner lobes, perhaps associated with 
a mixing layer. In addition, we have emphasised that simple com- 
pression by the bow shock should lead to enhanced RM's around 
the leading edges, but that the observed patterns depend on the pre- 
shock field. 

MHD simulations should be able to address the formation of 
anisotropic magnetic-field structures around radio lobes and to con- 
strain the initial conditions. In addition, our work raises a number 
of observational questions, including the following. 

(i) How common are anisotropic RM structures? Do they occur 
primarily in lobed radio galaxies with small axial ratios, consis- 
tent with jet-driven expansion into an unusually dense surrounding 
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medium? Is their frequency qualitatively consistent with the two- 
dimensional draped-field picture? 

(ii) Why do we see bands primarily in sources where the 
isotropic RM component has a flat power spectrum of low ampli- 
tude? Is this just because the bands can be obscured by large-scale 
fluctuations, or is there a causal connection? 

(iii) Are the RM bands suggested in tailed sources such as 
3C465 and Hydra A caused by a similar phenomenon (e.g. bulk 
flow of the IGM around the tails)? 

(iv) Is an asymmetry between approaching and receding lobes 
seen in the banded RM component? If so, what does that imply 
about the field structure? 

(v) How common are the regions of enhanced depolarization at 
the edges of radio lobes? How strong is the field and what is its 
structure? Is there evidence for the presence of a mixing layer? 

It should be possible to address all of these questions using a com- 
bination of observations with the new generation of synthesis arrays 
(EVLA, e-MERLIN and LOFAR all have wide-band polarimetric 
capabilities) and high-resolution X-ray imaging. 
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